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The physics of pions within a finite volume is explored using lattice regularized 
chiral perturbation theory. This regularization scheme permits a straightforward 

C3 ' 

computational approach to be used in place of analytical continuum techniques. 
CN . Using the pion mass, decay constant, form factor and charge radius as examples, it 

is shown how numerical results for volume dependences are obtained at the one-loop 

(N 

o 
o 



level from simple summations. 



^ i I. INTRODUCTION 

o 



Lattice QCD is one of the key tools for studying hadronic physics. It is a numerical tech- 
nique that employs a finite spatial volume, a finite extent in Euclidean time, and a nonzero 
spacing between sites on the spacetime lattice. Lattice QCD practitioners also choose un- 
physically large masses for up and down quarks due to the extreme cost of simulations at 
their physical values. 

The extrapolation to physical up and down quark masses can in principle be performed 
by using the low energy effective theory for continuum QCD, called chira. perturbation 
theory (ChPT).|2j The Lagrangian of ChPT contains an infinite number of terms, but to a 
specific order in the small chiral expansion parameters (for the pure pion theory these are 
m 2 I '(47T fir) 2 and p 2 /(47r/ 7r ) 2 with p being a small four-momentum) the number of terms is 
finite. ChPT has established itself as a valuable formalism for hadronic physics, and its use 
in connection to lattice QCD is just one important example. 

Similarly, the extrapolation in lattice spacing can be discussed within the effective theory 
for lattice QCD, which is simply ChPT extended to include the effects of the nonzero lattice 
spacing, a. This requires the addition of an infinite number of new terms to the continuum 
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ChPT Lagrangian, each of which is proportional to some positive power of a. To a specific 
order in the lattice spacing expansion, the number of a-dependent terms is finite and the nu- 
merical values of their coefficients can be determined in principle by matching to a particular 
definition of lattice QCD. Different lattice QCD Lagrangians (Wilson, Symanzik-improved, 
etc.) correspond to different lattice ChPT coefficients for the a-dependent counterterms. All 
of these additional terms become irrelevant in the continuum limit. 

Lattice ChPT can be defined within a continuum quantum field theory formalism, using 
(for example) dimensional regularization to handle ultraviolet divergences and retaining the 
lattice spacing only as prefactor for the a-dependent Langrangian counterterms mentioned 
above. Another option is to define lattice ChPT in an explicitly discrete spacetime.0, 
Q, The lattice spacing now plays the role of ultraviolet regulator in addition to being the 
expansion parameter for the a-dependent Lagrangian counterterms. With this approach, the 
lattice spacing appears explicitly in propagators and vertices and also in limits of integration 
for Feynman loop diagrams. The continuum and discrete methods are essentially equivalent 
when the inverse lattice spacing lies beyond the regime of ChPT (1/a > A x ~ m p ~ 4-7T f^) 
as is the case in typical lattice QCD simulations. One method or the other may be preferred 
for ease of use, or for theoretical discussions of the convergence properties of the ChPT 
expansion. [?], ls| 

The extrapolation in lattice volume within the framework of ChPT requires in general the 
inclusion of boundary- valued counterterms to the Lagrangian due to explicit boundary con- 
ditions except (as shown by Gasser and Leutwvler jol. \u\ ) for toroidal spacetime. In this case, 
the only effect of finite volume is the straightforward conversion of loop momentum integrals 
to loop momentum summations. For a review of recent finite volume ChPT calculations, 
see Ref. jl^. Some of the latest studies in the pion sector are those of Refs. 12, 13, 14]. 

In the present work, we explore the use of lattice regularized ChPT for computing volume 
dependences. The continuum limit must be identical to any viable continuum regulator, but 
lattice regularization has the feature of being easy to manage numerically. Beginning from 
a Lagrangian that displays the lattice spacing explicitly and also maintains exact chiral 
symmetry, [?| one can simply derive the Feynman propagators and vertices then type those 
directly into a computer program. Loop diagrams are just summations of a finite number of 
momentum values and the numerics are finite at every step. For a sufficiently small lattice 
spacing, observables must be independent of a. 



A brief preliminary discussion of this work can be found in Ref . Q| , but a more detailed 
study is presented below. Notation for the lattice regularized ChPT Lagrangian is estab- 
lished in Sec. |H] The computational method is introduced in Sec. IIHI by examining the 
two-point pion correlator. This gives the volume dependence of the pion mass, which repro- 
duces a result already known from continuum methods. 3, |l2| The two-point correlator also 
gives an explicit expression for wave function renormalization in the lattice regularized the- 
ory. Section IIVI contains a computation of volume effects on the pion decay constant which 
agrees with published continuum calculations. 0, Q New results are presented in Sec. |V| 
volume dependences of the pion form factor and the pion charge radius. Section IVT1 mentions 
some of the challenges that remain to be addressed if lattice regularized ChPT is to be em- 
ployed for the determination of volume dependences beyond the one-loop level. Appendix 1X1 
provides an explicit example of calculating Feynman rules from the lattice ChPT action, 
and Appendix [B] demonstrates the exact analytic agreement between volume dependences 
in dimensional regularization and in the continuum limit of lattice regularization. 



II. A DISCRETIZED SU(2) CHIRAL LAGRANGIAN 

The Lagrangian to be used in this work is an SU(2) version of the SU(3) meson Lagrangian 
introduced in Ref. Q|. Although only a few terms are presently required, here is the complete 
Lagrangian: 

C = C 2 + £ 4 , (1) 
C 2 = £ (vWtftvWtf) - £ ( x tc/ + xU i) , (2) 

~(h + k) {x ] U + X C/t) 2 + ^ 4 (V(%t v (%) (^U + x tf) - l 5 {FpflFfcj) 
~\h (F^U^U + F«V^UVtW) + l/ r ( x tU - xU^) 2 
~\{hi +h- h) (xh) + Uhu + fj (F^F% + F%F*) 

~(h -h 3 - k) (J^U + X U ] ) 2 + (x f U - X U ] ) 2 - 2 (x ] U X ] U + U^xU^x)) , (3) 
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where (...) denotes a trace, and summations over repeated Lorentz indices \i and v are 
understood. % is essentially the quark mass matrix, 

/ 



X = 2B 



(4) 



m u 
\ m d 

Throughout this work, we restrict ourselves to the isospin limit m u = = m q . We also 
choose the exponential representation for pions, 



U (x) = exp 



/ 



where r a is a Pauli matrix. The external fields are 

L^x) = exp [-iai^x)} = exp \-ia(V^(x) - A^x))} 
R p {x) = exp [-zar M (x)] = exp \-ia{V^(x) + A^x))] 

and the corresponding field strength tensors are discretized as follows: 



(5) 



(6) 
(7) 



Ma 2 F x 



ir „ 4 - X ll (x)X 1/ (x + a^Xfa + a v )Xl(x) 

-X l/ (x)Xj 1 (x - a M + a u )Xl(x - a M )X M (x - a M ) 

-X^(x - ajXKx - a M - a„)X M (a; - a M - a u )X u (x - a v ) 

- a I/ )X M (a; - a„)X„(:r + a M - a„)X^(x) (8) 

where X = L,R. As discussed in Ref. 0, a convenient way to avoid unphysical poles in the 
spectrum while maintaining invariance under parity is to use a nearest-neighbour derivative 
in the leading order Lagrangian, 

1 r- 



R^x)U(x + a^)LUx) - U(x) 



and a symmetrized derivative at next-to-leading order, 
1 



x 



2a 



R^(x)U(x + ajLUx) - RUx - a^)U(x - a^)L^(x - 



(9) 



(10) 



Notice that the Lagrangian in Eqs. (002) contains exactly the same number of terms as 
the continuum SU(2) ChPT Lagrangian[2j. As discussed in Sec. HJ the most general ChPT 
Lagrangian would contain extra terms proportional to positive powers of the lattice spacing. 
Since we are presently interested in volume dependences at the continuum limit, these extra 
terms are irrelevant and hence omitted for simplicity. 
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FIG. 1: Feynman diagrams contributing to the pion two-point correlator at one-loop level in ChPT. 



To conclude this section, we recall that the ChPT action is 



S = a 4 ££(*)-iWln 

T 1 ff \ 



2(1 - cos$(x)) 
$ 2 (x) 



where the second term is due to the integration measure, pj For SU(2) 



-2i 



f 



( n 7r 3 (a;) -tt 2 (x) \ 
n\x) 

y 7T 2 (x) — -n" 1 ^) y 



-7r 3 (x) 



(11) 



(12) 



III. THE PION MASS AND WAVE FUNCTION RENORMALIZATION 



The Feynman diagrams for the one-loop pion two-point correlator are shown in Fig. 
To evaluate them within lattice regularization, we choose a hyper-rectangular lattice with 
lattice spacing a in all four spacetime directions. The lattice is chosen to have N s sites in 
each of the spatial directions and N t sites in the temporal direction. Our goal is to consider 
the dependence of observables on spatial volume in the double limit a —>■ 0, aN t — ► oo with 
aN s held fixed. 

With Feynman vertices obtained from the Lagrangian of Sec. |H](see Appendix [X] for the 
derivations), the three diagrams of Fig. Q respectively become 



NLO 



~ x l ~ — XX 1 ~ cos ap^), 



2 n i ; \ 2/4X 2 ^ . 2 



(13) 
(14) 



F (6) = 1 

NLO WW^f 2 



^(ll2 + ha 2 x\ - 20 cos a Pv ~ 20 X cos afc P 



where 



+ 12 ^2 cos ap^ cos ak^ )D(k) 



D(k) 



a 2 x\ + 2 Z)u(l - cos afc M ) 



(15) 



(16) 
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is the pion propagator and 

x n = ^2Bm q (17) 

is the lowest-order pion mass in the continuum limit. The symbol "XV i n Eq. (fT3j) represents 
a sum over available lattice 4-momenta; for any function F, this means 

? fc2 ' h > h) = £ £ £ bv ^v7' ^v7' ^vT J • (18) 

Notice that the middle diagram in Fig. includes the measure contribution as well as the 
tree-level £4 contributions. 

The pion mass is defined as the energy of a stationary pion. The corresponding expression 
for the pion mass is the value of ip^ which solves T = when p — 0, where F is the sum of 
the three diagrams 

r = r L0 + rS + ra . (19) 

The result is 

M w = -arcsinh f^y 1 ) , (20) 

xi = 4 + 2 4h + 4 E ( \~ N *N7p D{k) + ° (a) ' (21 > 

J s t J 

Given numerical values for the Lagrangian parameters /, Bm q and Z3, the pion mass can now 
be computed directly from Eqs. pU|) and (}2"T]) for any lattice spacing and volume. As a — > 
the loop diagram diverges and these divergences are cancelled by the a dependence of the bare 
Lagrangian parameters /, Bm q and I3. (Since dimensional regularization retains no power 
divergences, / and Bm q would be scale invariant in that scheme. Lattice regularization 
does retain power divergences as a — * so the parameters / and Bm q do have a dependence 
in this scheme.) For any a ^ the loop diagram is finite, and for sufficiently small a the 
renormalized pion mass is independent of lattice spacing. 

To extract the volume dependence of the pion mass, one needs only the difference of M n 
at two different spatial volumes. The first two diagrams in Fig. ^ cancel in this difference 
leaving only the loop diagram. As a — > 0, the difference between two volumes must be finite 
because the only available Lagrangian counterterms were in the first two diagrams. The 
quantities and / appearing in the loop diagram are the leading chiral-order expressions 
for the mass and decay constant in the continuum limit. Following Ref. [l2j |. we employ 
/ = 86.2 MeV. One would expect results to become independent of lattice spacing for 
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FIG. 2: Fractional change in the pion mass as a function of spatial volume. 



a%l/ (47rf n ) ~ 0.2 fm, and we will choose N t 3> N s so that the temporal direction will not 
affect our extraction of spatial volume effects in any significant way. 

Figure El displays numerical results for the fractional change in the pion mass as a function 
of spatial volume, relative to the infinite volume pion mass, for M vr (L 00 ) = 100, 140 and 
300 MeV, corresponding to x n = 100, 142 and 321 MeV respectively. The computation at 
"infinite" volume, L^, is performed numerically simply by choosing a volume large enough 
to offer negligible deviations if the volume is increased yet further. Figure El shows explicitly 
the dependence of numerical results on changes to a, L^, and N t . As expected, heavier pions 
have an increased sensitivity to lattice spacing because loop integrals depend on the product 
ax n . The computation at a = 0.1 fm, = 8 fm and N = 5N S produces a fractional volume 
dependence for the pion mass that agrees with the known continuum result 0, ^ to within 
the resolution of this plot for the full range shown, 2 fm < L < 4 fm. As a confirming 
cross-check, this known continuum result is derived analytically from our lattice regularized 
expression in Appendix iBl 

In addition to the pion mass, Eq. (|19J) also leads to an expression for the wave function 
renormalization factor that will be required for all of the observables to be addressed below. 
Up to irrelevant lattice spacing effects, the two-point correlator can be parametrized as 



P 2 + xi + n 



-P 2 ) 



(22) 



r = 



Z + 




FIG. 3: Feynman diagrams contributing to the pion decay constant at one-loop level in ChPT. A 
wavy line denotes an axial vector current insertion. 



z2\ 



? + xi + EpO + {-v 1 - X 2 )Z'(X 2 ) + 5n-p 



where 



p 2 = 7^E si 



sin 



(23) 



(24) 



and <5£(— p 2 ) vanishes at least as quickly as (— p 2 — X 2 ) 2 for — p 2 — > X 2 . The renormalization 
factor, Z, is defined by 

".-.2 V2 1 

(25) 



r 



p 2 + x 2 



+ SE(-p 2 ) 



since X 2 = x 2 + S(X^) which leads to 

Z 



1 



1 - E'(X 2 ) 

One can read £(— p 2 ) directly from Eqs. (j 1311 15)1 by choosing p = 0, and this gives 

1 



/ 2 3NfN t a*f' 



^(5-3cosafc 4 ) J D(fc) +0(c 



(26) 



(27) 



IV. THE PION DECAY CONSTANT 



The three Feynman diagrams of Fig. |3] represent the three contributions to the pion decay 
constant up to one-loop order. Using the vertices and propagators from the Lagrangian in 
Eqs. (001), one finds the following expressions for those three diagrams, 



G L0 VZ = ^1 
a 

i2\[2 



sin apu + 2% sin' 



(?) 



G 



NLO 

(6) 

NLO 



(/ f 4^4 exp (^y^J cos sin(ap M ), 



(28) 
(29) 



—i 



sin ap M + 2% sin' 



faff/A 



v/2 



3W/ E(7-3cosa M D W . (30) 



V 2 / 

Choosing a stationary pion (p = 0) and inserting Eq. ()27jl for the wave function renormal- 
ization factor leads to 



<-*LO + ^NLO + ^NLO 



i^2 



f n [sin ap4 + 0(a)] 



(31) 
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FIG. 4: Fractional change in the pion decay constant as a function of spatial volume. 



where the one-loop pion decay constant is 



f ,, y , Xii2 - f T.(3-co S ak i )D(k) + 0(a) 

J s ^ J 



(32) 



In the difference between f n computed from two different lattice volumes, the first two 
terms in Eq. (|32|) subtract away. To this chiral order, the remaining parameters x w and / can 
be set to the (infinite volume) physical pion mass and decay constant. The resulting volume 
dependence of the pion decay constant is displayed in Fig. |3J The magnitude of the volume 
dependence is similar to that of the pion mass plotted in Fig. El but the sign differs — the 
decay constant is reduced as the volume shrinks, whereas the mass grows with shrinking 
volume. The computation at a = 0.1 fm, = 8 fm and N t = 5N S produces a fractional 
volume dependence for the pion decay constant that agrees with the known continuum 
result 0, Q| to within the resolution of this plot for the full range shown, 2 fm < L < 4 fm. 

The careful reader will notice that Fig. 0] has a slightly different normalization from the 
corresponding plot in Ref. 15]; this difference is higher order in the ChPT expansion, and 
is due to use of the physical mass and decay constant in Ref. [3] in place of the lowest- 
order parameters x n and /. Figures El and 0] of the present work show the familiar ratio, 
- ll/tM^/M^Loo) - 1] = -4, expected from Ref. fl. 
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FIG. 5: Feynman diagrams contributing to the pion electromagnetic form factor at one-loop level 
in ChPT. A wavy line denotes a photon. 

V. THE PION FORM FACTOR AND CHARGE RADIUS 



The pion electromagnetic form factor is obtained from the Feynman diagrams of Fig. 
and the charge radius can be extracted from the slope of the form factor at vanishing photon 
4- momentum. Using the Lagrangian of Eqs. (00)), the four diagrams evaluate as follows, 



Hi 



H 



LO 



(a) 
NLO 



2Z (-iaqu\ . ( p + p' 
— exp [ — ) sin a 



a ~ \ 2 ) \ 2 
2kxl 



(33) 



exp — — - cos — - sin a{p + p )[_ 



H 



H 



(6) 
NLO 

(c) 
NLO 



a 3 / 2 
-10 



exp 



af 2 ~' r V 2 ) V 2 ) 

2l R f-iaq^\ ^ r . . , . , . 

— - I [sin ap^ sin ap u — sin ap^ sin ap v 



sm cos I -y 1 j ,(34) 



3N^N t a 3 f 
4 



exp 



sin 



(35) 



y , V (/ i/2 exp (^^J I>na(£: + g/2) M £ cosa(p - k) u D(k)D(k + q), (36) 

where p is the incoming pion momentum, q and /J are the incoming momentum and Lorentz 
index of the external photon, and p' = p + q. The contribution from H^l can be simplified 
by removing terms that are odd under interchange of k and — (k + q), since these vanish 
after summation over k. The result is 



(c) 
NLO 



exp 



^sina^^j Y, sin2 a{k + q/2)^D{k)D{k + q). (37) 



N^N t a 3 f 2 ~' r \ 2 J V 2 
The pion form factor, F(q 2 ), can be obtained explicitly by choosing /x = 4 as follows 



it , rr(a) , rr(6) , tt(c) 
-"LO + -"NLO + -"NLO + -"NLO 



-F(q)exp — ■■- 
a V 2 



sm a 



(38) 



which gives 



i^ 2 



,2 f 2 



[sin ap 4 sin ap^ — sin ap^ sin ap u ] sin ag^ cos ( 
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FIG. 6: Fractional change in the pion form factor at Q 2 = M 2 as a function of spatial volume. 



N*N t a 2 f 2 ^ 
+0(a). 



^2 cos ak^D^k) + 



(39) 



It is interesting to consider the q — > limit, since vector current conservation should require 
F(0) = 1. From Eq. (J39|) . we see that the term containing Z 6 does vanish in the q — > limit. 
The cancellation at g = of the two summation terms from Eq. (J39|) is easily demonstrated 
in the notation of a temporally infinite lattice, 



^3 E f° 7^ (2 sin 2 ak A D 2 (k) - cosak 4 D(k) 

sin akiD(k)) 



1 /*/« dfc 4 d 

0. 



(40) 



Fig. El shows the numerical results for the volume dependence of the pion form factor at 
q 2 = M 2 (meant to represent a typical ChPT mass scale) obtained from Eq. (|39|) . where for 
numerical ease we work in the Breit frame. As is evident from the plot, the form factor's 
fractional volume dependence has a similar magnitude to that obtained for the pion mass 
and decay constant. 

The pion charge radius is extracted from the slope of the form factor at q 2 = 0. Choosing 
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FIG. 7: Fractional change in the pion charge radius as a function of spatial volume. 

q = (0, 0, q 3 , 0) for deflniteness, we find 
(rl) = -6 1im^ 

= — + Ar , Ar f , y sin 2 aki cos ak 3 D 3 (k) — 4 sin 2 afc 4 sin 2 ak 3 D 4 (p) (41) 
/ N*N t f 2 Y 

A graph of the fractional volume dependence of this quantity is provided in Fig. The 
magnitude of the effect is dramatically larger than for the mass, decay constant and form 
factor simply because the charge radius has (volume-dependent) loop contributions at its 
first nonzero ChPT order. 



VI. SUMMARY AND OUTLOOK 

The pion mass, decay constant, form factor and charge radius have been computed from 
0(p A ) chiral perturbation theory in a finite volume by using lattice regularization. A sug- 
gested advantage of this regularization scheme is that the renormalization can be carried out 
numerically, leaving fewer analytical steps to be performed. Explicit expressions for these 
observables are given as four- dimensional finite sums in Eqs. (jZH), (JB2J), (jHEJ) and (jUJ). 

The dimensionally regularized expressions for thepion mass and decay constant are known 
to be one-dimensional sums over Bessel functions 2, [l3|, and results from the two regular- 
ization schemes agree numerically. In essence, dimensional regularization arrives at a more 
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compact result (i.e. fewer summations) because in that method more of the renormalization 
is done analytically. Given the small computational cost of the four-dimensional summa- 
tions, the lattice regularized result is also quite usable in practice. In addition, Appendix iBl 
demonstrates that the continuum limit of the lattice regularized result is analytically identi- 
cal to the dimensional regularized expression, if one chooses to complete the entire analytical 
calculation instead of the computational scheme proposed in this work. 

As emphasized in Ref. jll|, it is necessary to extend discussions of volume dependence 
to the two-loop level, and perhaps beyond, so that the rate of convergence can be explored. 
In general, two-loop renormalization is substantially more involved than the one-loop case 
so the reduction of analytical effort obtained by using lattice regularization could be of 
considerable value. The extension of lattice regularization to two loops will involve the 
determination of numerical values for the Lagrangian's low energy constants (and their scale 
dependences) since they will no longer subtract away in the difference between two volumes. 
It will also require an understanding of the interplay between power divergences, l/a n , 
and volume dependences, l/(aN s ) n . In particular, one does not want to rely on numerical 
cancellations among diverging summations. These issues are currently under investigation, 
in hopes of extending this practical computational method to the domain of multi-loop 
ChPT calculations. 



Acknowledgments 

The authors thank Daniel Mazur for his involvement in the initial stages of this research, 
Georg von Hippel for a careful reading of the manuscript, and Hermann Krebs for useful 
discussions. This work was supported in part by the Deutsche Forschungsgemeinschaft, the 
Natural Sciences and Engineering Research Council of Canada, and the Canada Research 
Chairs Program. 



APPENDIX A: LATTICE FEYNMAN RULES FOR THE PION MASS 

This appendix contains an explicit derivation of the Feynman rules required for a one- 
loop computation of the volume dependence for the pion mass. Only the third diagram of 
Fig. ^ contributes, which contains the pion propagator and the four pion vertex. 
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Beginning from Section |H] and using ellipses to denote terms that do not contribute to 
the pion two-point function, the relevant terms in the action are 



SS = a 4 ^£ 2 (x) + ... 



a 2 / 2 



X)Tr \u\x + fi)U{x) + U\x)U(x + fi) 



Bm q Tr \u(x) + U\x) 



X,fJ, 



+ 



a 4 x 2 



-a 2 6 ab J2 K a {x) [n\x + jl) - ir\x)} + ^5 ab £ 7r a (x)7r 6 (x) 



(Al) 



X,fl 



with /2 = ae^. This will now be expressed in terms of the Fourier transform, 



1 



71 (X) = 



Em 



Ak-x 



where the summation extends over the set of momenta 

/27rni 27rn2 2imz 

= V oivT' aiv7' aiv7' ~^/vT 



(A2) 



(A3) 



with 



rtj 



" 2 ' 2 + 2 ' 

— + 1 — — 1 

2 ' 2 ' "' 2 ' 



(A4) 
(A5) 



and we have chosen A^, N t to be even, in order to keep the presentation simple. Using the 
relation 

Jk-x _ AT 3 ]\t a( 4 ) 



leads to 



5S 



a 2 5° 



v , v E C^W) 



Ak'-fl 



a 4 xl5 ab 



4 raf> 



a 4 <5 



E 



' + ^E( 1 "CosA;- / i) 



^r a (fc)^(-£;) + . . . 



(A7) 



The Euclidean two-point correlator for incoming pion fields 7r a (j») and tt (q) is 



d 



d 



55 



_ra6^(4) 
°p+<?,0 



4 + — E( x - cos P'P) 



(A8) 



a \N$Nt ) d \N§Nt ) V iv « iv */ 

The pion propagator is the negative inverse of this expression with p = —q and a = b, and 
is a 2 D(k) of Eq. (HSJ). 
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To derive the four pion vertex, we again begin from Section |H] and use ellipses to denote 
terms that do not contribute, 



SS = a 4 ^£ 2 (x) + ... 



a 2 / 2 



Y^ \u\x + fi)U{x) + U\x)U(x + fi) 



Bm q Y Tr \u{x) + U\x) 



X41 



+ 



;Tr(rVVV) Y [* a ( x + A) 71-6 + + /i)vr d (a; + p) 



48/ : 

-4vr a (x + fl)ir b (x + /t)vr c (x + fl)iT d (x) + 6vr a (x + /x)tt 6 (x + fi)iT c (x)iT d (x) 
-47r a (x + /i)7r 6 (a;)7r c (x)7r d (x) + 7r a (x)7r 6 (2;)7r c (x)7r d (x)" 

;Tr(T a T b T c T d )Yn a (x)n b (x)TT c (x)ir d (x) + ... 



4^2 



ax 



48/ 2 



- , Rf2/vl2/v4 Tr(rVW) £ 7f a (A;)7r 6 (^)^ c (^)^ d (^)^ (A+fc ' +fc " +,s '"^ 

40/ 1V S l\ t x,fi,k,k',k",k"' 
e i(k+k'+k"+k"')-fi _ ^ e i(k+k'+k")-fjL _|_ Q e i{k+k')-fi _ ^gifc-A _|_ ^ 



a 4 x 2 



48/ 2 iV s 12 iV t 
.,2 



^Tr(T a rVr d ) ]T ^ a (A;)^(A; / )r(A; /, )^(A;'")e i(fc+fc ' +fc " +fc "' ) - :c 



r U Ul h" h !!! 



rTr(r a rVr d ) £ n a (k)n b (k')r(k")n d (-k — k' - k" 



k,k',k h 



2Af 2 NlN? 



(A9) 



The Euclidean four-point correlator (i.e. the Feynman rule) for incoming pion fields vr a (p), 
7r 6 (g), 7r c (r) and n d (—p — q — r) is 

d d d d ( 5S \ 



N?N t ) d \NlN t ) a \NlN t ) a { N*Nt ) V s * 



2 



3a 2 / 2 



(5 ab 5 cd + 5 ac 5 bd + 5 ad 5 bc ) Y (1 - c °s?V - cos g M - cos r M - cos(p + g + 



a 2 / 2 



+^5 ac o fed ^ cos(p + r)„ + -4^S ad 5 bc Y cos(g + 



' f- 1 ' 



(A10) 



The a — > limit reproduces the standard continuum Feynman rule as expected. 
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APPENDIX B: ANALYTIC DERIVATION OF THE CONTINUUM LIMIT 



In the main body of this article, lattice regularized results were left in the form of loop 
summations over products of Feynman rules, since this is sufficient to produce numerical 
results. As expected, the numerics agreed with analytic dimensional regularized calculations 
where available, since physics does not depend on regularization scheme. Using the volume 
dependence of the pion mass as an explicit example, this appendix verifies that continuing 
the analytic steps in the lattice regularization approach, and taking the continuum limit, 
leads to the same analytic result that comes from dimensional regularization. 

From Eqs. (|2*nj) and (|21j) . lattice regularization produces the following difference in pion 
mass for two lattice volumes, 



[M^L)-M n {L r ) 



I finite a 



lim 

Nt^oo 



E 



3 — 2 cos aki 



Aa 2 f 2 N^N t ^ a 2 x\ + 2 £ M (1 - cos ak^ 
x n ^ 3 — 2 cos ak A 



E 



Aa 2 f 2 (N' s fN t a 2 x% + 2 - cosak'J 
where the summations extend over the set of momenta 



+ 0(a), (Bl) 



k 



2-ktii 2irn2 2im% 27m 4 \ 
aN s ' aN s ' aN s ' aN t ) 



k' 



27m[ 27m' 2 2irn f 3 27rn 4 



with 



N 



rtj 



n ■ 



n>4 



N 

2 ' 2 

Nt n: 

~2 



- - + 1 

2 ' 2 + ' 

_iV £ _Nt 



aN' s ' aN' s 



2 

N 



aNL ' aN t 



(B2) 



N t 



- 1 



(B3) 
(B4) 
(B5) 



and N s ,N' s ,N t} Nj. even. 



Using a double angle formula from basic trigonometry, Eq. (jBl|) can be re-expressed as 

x n ^ 1 + 4sin 2 (afc 4 /2) 



[M n (L) - M^L')] 



finite a 



lim 

Nt^oo 



Aa 2 f 2 N^N t ^ a 2 x\ + 4^ sin 2 (a£; M /2) 
x n ^ 1 + 4sin 2 (a£; 4 /2) 



E 



O(a). (B6) 



Aa 2 f 2 {N' s fN t ^ a 2 x\ + AJ2^m 2 (ak^/2)_ 

Since this mass difference is finite even in the continuum limit, and since the Taylor expansion 
of sin# satisfies absolute convergence term by term over the entire range of interest, —tt/2 < 
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9 < tt/2, the leading volume dependence is obtained by retaining the leading term in this 
Taylor expansion. The result is 



[MM - M«(L')] mte a 



, f'9 

4a 4 / 2 Nt^oo 
+0(a). 



1 



£ 



i 



£ 



i 



(B7) 



As stated in Sec. IIII[ our goal is to consider the dependence of observables on spatial volume 
in the double limit a — > 0, T = aN t — > oo with L = aA s (and for now L' = aA 7 ^ also) held 
fixed. With this simple change of variables (but leaving a finite momentarily), we obtain 



[M W (L) - M^L')}^ a = lim 



£ 



4/2 T— »oo [L 3 T V 4 + E,(2vrn,/L) 2 + (27m 4 /T) 5 



.A3T 1 ^ 



(L') 3 T ^ xl + E ? (2vrn;./L') 2 + (27m 4 /T) 2 



0(a) 



X T 



— V 

T3 ^ 



4/ 2 y-^r/a 27T [I, 3 ^ £ 2 + Y,j(^ nn j/ L) 2 + £> 4 



1 



ka£ 



+ 0(a). 



(B8) 



(L') 3 V4 + E J (27r^./L') 2 +p|. 

We now take the continuum limit, which merely extends the bounds of summation and 
integration to ±00. The integral over p 4 can be performed in closed form, 



M^L) - M W (L') 



X -n 



]f 2 



£ 



1 



1 



A3 £ 



1 



L 3 V Jxl + E^Trn./Ly (L'f ^ l x 2 + E .( 27m ;./L') 



(is- 



L 3 



£< 



'•E,( 2 ^/ L ) 2 1 



A3 



^M,/L')> 



(B9) 



Letting L' — > 00 gives 
M W (L) — M n — 



x T , 



00 g - 
ds- 



]f 2 ^ Jo y/s 



f 2 ^Jo y/s 



1 

I 3 

1 

I 3 



£ e - 



d 3 p 



-sp 



£< 



;J2 A^nj/Lf _ 



(4tts) 3 / 2 _ 



(BIO) 



where M n = M n (oo). We can now make use of a relation that appears in Ref. 



• [16]. 



£ 



00 

71 — 7r 2 n 2 /r 



£ 



(Bll) 
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to obtain 



M„(L) - M n 



X 7] 



64/ 2 7T 2 

x n 
64/ 2 vr 2 

16/ 2 7T 2 L 



ds- 



J2 e -n 2 L 2 /(4s) _ 1 



oo g" 

ds — 



-n 2 L 2 /{As) 



E 



K^x^LJnl + n| + nl) 



n 2 + 7^2 + «| 



(B12) 



where ifi(x) is a Bessel function of the second kind. This is the result known from dimen- 
sional regularized calculations. Though expressed triple summation over ni, 712 
and ri3, the function only contains the sum of squares, n^+n^+nl, thus allowing M n (L) — ML 
to be represented by a one- dimensional summation when multiplicity factors are defined. [12( 
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